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ABSTRACT 

We present a study of the molecular CO gas and mid/far infrared radiation arising from the 
environment surrounding the Wolf-Rayet (W-R) star 130. We use the multi-wavelength data 
to analyze the properties of the dense gas and dust, and its possible spatial correlation with 
that of Young Stellar Objects (YSOs). We use CO J=l-0 data from the FCRAO survey as 
tracer of the molecular gas, and mid/far infrared data from the recent WISE and Herschel 
space surveys to study the dust continuum radiation and to identify a population of associated 
candidate YSOs. The spatial distribution of the molecular gas shows a ring-like structure very 
similar to that observed in the H i gas, and over the same velocity interval. The relative spatial 
distribution of the H i and CO components is consistent with a photo-dissociation region. 
We have identified and characterized four main and distinct molecular clouds that create this 
structure. Cold dust is coincident with the dense gas shown in the CO measurements. We 
have found several cYSOs that lie along the regions with the highest gas column density, and 
suggest that they are spatially correlated with the shell. These are indicative of regions of star 
formation induced by the strong wind and ionization of the WR star. 

Key words: Stars: individual: WR 130 - ISM: bubbles - Hii regions - Infrared: ISM - Stars: 
formation 


1 INTRODUCTION 


Massive 0-type stars and their descendants have a major impact on 
their surrounding interstellar medium (ISM) via their high through¬ 
put of ionizing photons, their energetic winds, their expanding H ii 
regions, and their final explosion as supernova. 

During their evolutionary courses, massive stars possess quite 
different mass loss rates and wind speeds, generating a vari¬ 
ety of str uctures such as bubbles, shells (or multi-shells) in 
the ISM (Castor. McCray & Weavei] 1 19751 : IWeaver et~^ Il977l : 
iGarcia-Segura & Mac Lowlll995h . Due to the high number of ul¬ 


traviolet photons emitted by the star during the 0-type p hase, an 
Hii region is formed around the star dWeaver et al.ll 19771) . As this 
hot region expands, neutral material accumulates between the ion- 
izati on front and the shock fr ont that precedes it in the neutral gas 
fsee lDvson & Williamsil 19971) . Simultaneously with the ionization 
of its surroundings, stellar winds are blown by massive stars, which 
may affect the evolution of the H ii region. 

The following phase in the ev olution of the star depends on 
its initial mass fe.g. ICrowthei]l200% . For masses higher than ~ 60 


Member of the Carrera del Investigador Cientifico of CONICET, Ar¬ 
gentina. 

t Post-Doc Fellow of CONICET, Argentina 


M© the star could go directly to a Wolf-Rayet (W R) phase, with- 

a nng first throught a RSG or LBV phase (e.g. ISmith & Contil 


During the WR phase, the stellar wind strengthens (could be 
an order of magnitude higher than in the main sequence phase) and 
may clear the region around the ionizing star, creating a stellar- 
wind bubble filled with high temperature (T ~ 10^“^ K) and low 
density (n ~ 10"^ - 10"^ cm“^) shocked gas. In due time, this hot 
cavity may drive expanding gas shells that may eventually inter¬ 
act with the shell related to the expanding H ii region. In this way, 
the supersonic winds can entrain and accelerate the surrounding 
gas, thereby injecting momentum and energy into the surrounding 
gaseous environment, strongly affecting its evolution. 


The structures created by the action of massive stars, usually 
striking at radio frequencies, are also very bright at far infrared 
wavelengths, due to the presence of dust grains that absorb the 
UV stellar photons and cool down emitting in the infrared. The so- 
called infrared dust bu bbles can now been studied i n detail since the 
realease of the WISE jWright et al.ll201^. Spitzer (iBeniamin et al.l 

(iMolinari et al.ll2 


(e.g. 


12003: ICarev et al.|2009h. and Herschel (iMolinari et al.ll2010h data 


Zavagno et al .1120 1 ol: IPaladini et al .1120 1 2h . 


One important consequence of the interaction between masive 
stars and their local ISM is that the physical conditions present in 
the dense expanding shells are ideal for star formation activity. Sev- 
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eral studies in fact show the presence of young stellar object can¬ 
didates (cYSOs) wh ose birth was probabl y triggered by the expan¬ 
sion of Hii regions (IZavagno et al.ll201Ql. and references there in). 
However, similar studies towards Wolf-Rayet stars, r epresentative 
of more evolved stages of massive stars, are scanty dCappa et al.l 
l2010l:lLiu et allllOialLiu. Wu & Zhangll2013h . Thus, it is still un- 
clear whether these older bubbles can create regions of active star 
formation around them. 

In this work we present the second part of a multiwave- 
length study of the surroundings of WR130, a WN8(h) s tar 
(Ivan der Huchtll200ih or, as proposed bv ISmith & Contil (l2008h . a 
WNH. The place th at this type of WR stars, w hich are the least 
chemically evolved JConti. Leep & Perrvlll983h . has in the evolu¬ 
tion sequence of a massive star is not completely understood. Based 
on the mass discrepancy observed between WNH and H-poor WR 
stars (the WNH stars have higher masses), and on the fact that 
the most luminous W NH stars tend to be more H-rich than LB Vs, 
ISmith & Contil (l2008h concluded that the WNH stars are pre-LBVs, 
and that via the violent eruptions of the LBV phase the star removes 
the hydrogen envelope and becomes a WN star. In this scenario, 
WR 130 could be in a pre-LBV phase. The observed impact that 
WR 130 has on its environment seems to support this possibility. 
Indeed, a very well defined ring-like Hii region, first noticed by 
iHeckathorn, Bruhweiler & Gullf(ll982h . is associated with this star. 
Later, ICichowolski et iD|2^ studied the region using a multi¬ 
wavelength approach that involved DRAG H i line data and contin¬ 
uum at 408 and 1420 MHz, HllOof VLA data, HIRES 60 and 100 


jjim data dPowler & Aumannlll994 ), and 12 an d 25 jjim Mid-Infrared 
Galaxy Atlas images dKerton & Martinl200Cll) . From all the data an¬ 
alyzed, [Ochowo^^e^lJiliQoIFobserveda well defined ring-like 
structure, which is consistent with a scenario were the action of the 
O star was directly followed by the action of the WR star, since the 
material expected to be present if the LBV phase took place is not 
detected. 

In this paper we now present a study of the CO clouds detected 
in the area of WR 130 and analyze the properties of the dust making 
use of the new available WISE and Herschel data, which allow us 
to provide new insights into the gas and dust properties of the inter¬ 
stellar medium surrounding a Wolf-Rayet star. We also search for 
the presence of triggered star formation towards the surroundings 


ofWR130. 


2 OBSERVATIONS 

In this paper, we use the 1420 MHz radio continuum image of the 
WR 130 area from the Canadian Galactic Plane Survey (CGPS - 
iTavlor et ^ l2003h . Because this image is the result of a mosaic 
containing more than one pointing field, it has a slightly better 
signal-to-noise ratio and dynamic range. It also has better flux and 
position registration, and a slightly better resolution, factors which 
are important when estimating the fiuxes of some discrete sources 
seen projected over an inhomogeneous background. 

At infrared wavelengths, we have used the data r etrieved from 
the W ide-Field Infrared Survey Explorer (WISE) dWright et alJ 
bold) , which is an all-sky survey that contains images in four in¬ 
frared channels centred at 3.4, 4.6, 12, and 22 fim. These channels, 
respectively, have angular resolutions of 6.1, 6.4,6.5, and 12 arcsec. 
We obtained all the images from the NASA/TPAC Infrared Science 
Archive (http://irsa.ipac.caltech.edu). __ 

We have also included Herschel data dPilbratt et al.l l201(]h ) 
from the Herschel Infrared Galactic Plane Survey (HiGAL, 


iMolinari et aljboioh obt ained in its five phot ometric bands at 70 
and 160 fim with PA CS dPoglitsch et ^bOld) . plus 250, 350, and 
500 jjm with SPIRE dGriffin et al.ll2010l) . The maps have been pro¬ 
duced by the improved version of the ROMAGAL pipeline us¬ 
ing UNIMAP, and therefore, are optimized to preserve the ex¬ 
tended/diffuse e mission, plus to avoid artifacts arou^^nd br ight point 
sources (see e.g. lTraficante et alJl201ll : l^azzo et alJ2015h . The Hi- 
GAL maps have angular beams with a FWHM of about 5, 12, 18, 
25 & 36" at 70, 160, 250, 350 & 500 //m, respectively. And the 
maps themselves were created with pixel sizes of 3.2, 3.2, 4.5, 6.0, 
8.0, and 11.5" at 70, 160, 250, 350 & 500yum, respectively. 

Molecular CO / = 1 - 0 observations were obtained from the 
Five College Radio Astr onomical Observa tory (FCRAO) CO Sur¬ 
vey of the Outer Galaxy dHever et al.ll998h . The angular resolution 
is 45" , the velocity resolution 0.127 km s ~^, and the rms noise 0.2 
K 


3 CO EMISSION 

Figure[T]shows a set of images of the CO(l-O) emission distribution 
within the velocity range from about -11 to 3 km s~^ (all velocities 
are with respect to the Local Standard of Rest (LSR)). Each image 
is an average of the emission in ten consecutive velocity channels, 
where the corresponding central velocities are indicated in the top 
side of each box. The velocity range show n was chosen accord¬ 
ing to the work of ICichowolski et ^ (l200ll) . where they found an 
Hi structure related to WR 130 in the velocity range from -10.3 
to 1.2 kms~^ To indicate t he location of the Hii region found 
by ICichowolski et al.l (1200ll) . a 1420 MHz contour delineating it 
is shown in the top left panel. 

As can be seen in Fig. [T] several molecular clouds are present 
in the area of WR 130. In particular, at -4.65 km s“^ the star lies, 
in projection, inside a well defined molecular shell. To compare 
the molecular di stribution with the ioni zed and atomic gas re¬ 
lated to WR 130 dCichowolski et al.ir200ll) . in Fig.|2]the CO emis¬ 
sion averaged between -11 and 3 kms“^ is shown, along with 
the 1420 MHz and Hi emissions. The left panel of the figure 
shows the emission distribution at 1420 MHz. The radio contin- 
uum emission in the area of WR 130 was extensively analyzed by 
ICichowolski et al.l (l200ll) . who concluded that the ring like feature, 
named G68.1.1-hl.l, is the radio counterpart of the optical region 
Sh2- 98 and it is related to WR 130. On the contrary, based on RRL 
data, ICichowolski et al.l (1200ll) found that the brightest compact 
source, named G68.14-^0.92, observed at (/, Z?) = (68U4, 0?92) is 
an Hii region located at a distance of about 12 kpc, and thus unre¬ 
lated to WR 130. 

In the middle and right panels of Fig.|2l the CO and H i emis¬ 
sions are shown, respectively, with the radio continuum level de¬ 
lineating the Hii region superimposed. From these images it is 
clear that the ionized gas is surrounded by an atomic shell and 
that several molecular clouds, which have probably formed part 
of the molecular cloud where WR 130 was born, are still present. 
The observed emission distributions clearly indicate that a photo- 
dissociated region (PDR) was formed around WR 130. 

To estimate the molecular mass of the clouds probably related 
to WR 130 we first estimate the column density through the empir¬ 
ical relation 


N(H2) = XxIco 

where X is the CO-to-H 2 conversion factor, for which we consid- 
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Figure 1. CO emission distribution in the velocity range from -11 to 3 kms“L Each image is an average of ten consecutive velocity channels, where the 
corresponding central velocities are indicated in the top side of each box. Contour levels are at 0.3, 0.7 and 1.1 K. Angular resolution is 2L Beam shape is 
shown in the upper left corner of the bottom left p anel. The 9.5 K level at 1420 MHz is indicated in the top left panel to delineate the ionized region found to 
be related to WR 130 bv ICichowolski et ahl (120011) . The star symbol indicates the location of WR 130. 
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Figure 2. Left panel: CGPS 1420 MHz emission distribution as shown in ICichowolski et aP j2onih . Contour level is at 9.5 K. Middle panel: CO emission 
averaged between -11 and 3 kms“^ Contours levels are at 0.3, 0.7 and 1.1 K. The dotted contour corresp onds to the 9.5 K emission at 1420 MHz. Right panel: 
CGPS Hi emission distribution averaged between -10.3 and 1.2 kms“Has in ICichowolski et all (1200 ih ) with the 1420 MHz and CO contours superimposed. 
In all the panels the star symbol indicates the location of WR130. 


ered the varying relation w ith the Galactocentric radius R given by 
iNakanishi & Sofuel (l2006h . 

X[cm-^ K-^ kms-^] = 1.4 x lO^^ exp(R/ll kpc), 

and Ico is the integrated line emission in the velocity range where 
the cloud is observed, Ico - jTdv. 

Figure |2] shows a comparison of the emission observed in 
the radio continuum at 1420 MHz, the CO and Hi emission, the 
latter two being averages over a few velocity channels. Assum¬ 
ing that the molecular gas show n in Fig. |2] (middle panel ) is re¬ 
lated to the HI feature found by ICichowolski et al.l (l200lh as be¬ 
ing created by WR 130, we adopt for the m olecular clouds the 
distance inferred bv ICichowolski et ^ (1200 ih . i.e. D = 5 kpc. We 
thus obtain that the structure is located at R = 8 kpc, which yields 
X = 2.9 X 10^^ cm"^ km s“^ Then, we estimate the total molec¬ 
ular mass related to WR 130 based on the averaged emission shown 
in Fig. 121 as well as the mass of the four individual molecular clouds 
indicated in Fig. [T] These clouds were selected because they host 
infrared source candidates to be young stellar objects, as will be 
discussed in Section |5] The molecular masses were derived from 
Mq = 4.2 X N(H 2 ) A, where D is the distance in pc and 
A is the solid angle in steradians. Assuming an elliptical geometry 
for the individual clouds , the volume density can be estimated as 
n(H 2 ) (cm"^) =10 Mo/a where a and b are the major and minor 
semiaxes, respectively, in pc. The masses and densities obtained, 
along with the velocity intervals considered for each estimate are 
shown in Table [T] 

We have also estimated the dust masses associated with each 
cloud, which are given in the last column of Table [T] They were 
obtained using the Herschel data. For each cloud fluxes were mea¬ 
sured at 70, 160, 250, 350, and 500 jim by rebinning all the images 
to the CO resolution and subtracting emission from the local back¬ 
ground. Then the spectral energy distribution (SED) of each region 
was fitted using the models created by DUSTEM (see Section lU 
allowing the estimation of the dust masses associated with each re¬ 
gion. 

The kinetic energy stored in the CO shell can be estimated 
as Ekin = 0.5 Msheii Fgxp, where Vexp is the expansion velocity of 
the shell and Msheii is the total (molecular, atomic and ionized) 
shell mass. Adopting an expansion velocity equal to half the ve¬ 


locity interval where the structure is observed, Vexp = 7.0 +1.3 
km s"^, the molecular mass given in Table [T] and the atomic and 
ionized masses estimated bv Tcichowolski et^ (l200lh . 1450 and 
3000 Mo, respectively, we obtain = (2.5 + 1.0) xlO"^^ erg, as¬ 
sumin g a 40% error for the masses. Although ICichowolski et al.l 
(l200lh concluded that WR 130 could have alone created the ob¬ 
served structure, it is important to note that they did not take 
into account the molecular mass present in the shell, which con¬ 
siderably increases the kinetic shell energy. Thus, we can com¬ 
pare now the new value obtained for with the mechanical 
energy deposited in the ISM by the wind of the Wol f-Rayet star, 
E^J^ - (0.7 - 2.2) X 10^^ erg (ICichowolski et aklboOlh . We obtain 
6 = Ekin/Em = 0.007 - 0.5. The ratio e measures the energy conver¬ 
sion effic iency in the shell, and according to evolutionary models 
6 < 0.2 (IKoo & McKeelll99^ . Thus, not all the possible values 
of 6 are compatible with the scenario where the energy injected 
during the Wolf-Rayet phase is enough to create the structure. In 
this case, the contribution of the energy injected during the O star 
phase and/or other massive stars, should be considered. As men- 
ti oned in the Introductio n, WR 130 is a WNH star, and according 
to ISmith & Contil (l2008h its age would be of about 2-3 Myr and 
its initial mass of at least 60 M©. A rough estimation of the en¬ 
ergy injected by such a star during its main sequence yields Eu, - 
(2.5 - 3.5) X 10^^ erg dde Jager. Nieuwenhuiizen & van der HuchtI 
Il988h . which would be enough to create the observed structure. 
We have nevertheless looked for the presence of other massive 
stars in the region. We qu eried the available catalogue s such as the 
Galactic 0-Star Catalog dMaiz Anellaniz et ^ l2013h. the Early- 
Type Emission-Line Stars Catalogue (Iwackerlindll970h . the Cat¬ 
alogue of Be stars dJaschek & Egret Il982h. the H-alpha Stars i n 
the Northern Milky Wav Catalogue dKohoutek & WehmeveJ 19971) . 
and the Catalog of Galactic OB Stars dReedI12003 ). for early-type 
and emission stars. No stars were found in any c atalogue. The only 
massiv e star located near by is, as mentioned bv ICichowolski et aP 
d200lh . an OB st ar, which has an uncertain spectral type and no 
distance estimate dStock. Nassau & SteDhensonlll960h . It is located 
in projection not in the center of the st ructure but onto the shel l 
(there is a second OB star mentioned bv ICichowolski et ^ d200lh 
but its location is actual ly outside the structure, see Eigure 1 of 
ICichowolski et al.l d200lh ). Although we can not completely rule 
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Table 1. Physical parameters of the molecular gas. 



01 

(kms“^) 

V2 

(kms“^) 

Mass 

(lO^Mo) 

Density 
(10^ cm“^) 

Dust Mass 

Mo 

Total 

-11 

3 

46 

- 

- 

Cloud 1 

-3.5 

1.84 

2.2 

1.2 

114 

Cloud 2 

-2.35 

4.25 

4.6 

2.4 

180 

Cloud 3 

-10.5 

-5.0 

2.6 

0.8 

125 

Cloud 4 

-7.18 

-0.7 

6.0 

0.8 

252 



Galactic Longitude 

Figure 3. WISE 12 jum image of the region. The contour levels correspond 
to the averaged CO emission between -11 and 3 km s“^ and are at 0.3, 0.6, 
and 0.9 K. 

out the possibility that the OB star may be playing a role in creating 
the shell structure, we think that the action of WR 130 is sufficient 
and most likely dominant in the region. 


4 INFRARED EMISSION AND DUST PROPERTIES 

In this section we make use of the WISE and Herschel data to ana¬ 
lyze the distribution of the dust as well as its properties. 

As shown in Section [3] the ionized region related to WR 130 
is interacting with the molecular gas. As a consequence a PDR is 
created in the interior layers of the molecular clouds, which can be 
seen as polycyclic aromatic hydrocarbon (PAH) emission in WISE 
A at 12 jLim. Figure [3] shows that the 12 jam emission displays a 
cavity surrounded by a roughly thin annular shell emission. 

Figure |4] shows the distribution of the warm and cold dust 
traced by the 12 jum W3 band (PAHs; blue), 22 jdm W4 band 
(VSGs; green) and Herschel 250 jam (Big Grains; red). Overall the 
spatial distribution of the three components looks quite similar, at 
least at the angular resolution of the WISE data, 6 and 12 arcsec- 
onds for W3 and W4, respectively. The biggest difference takes 
place on the rim near the WR 130 source itself, where a warmer 
dust component appears closer to the source, suggesting a more 
typical photodisociation structure, with the warm gas facing the 
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Figure 4. Color composite image of the area of WR 130. Blue color shows 
the emission at 12 jum (W3 WISE band), green color represents the emission 
at 22 jum (W4 WISE band), and red color shows the emission at 250 jum 
(SPIRE). The location of WR 130 is indicated by the white circle. 

ionizing source. WR 130 itself, is detected by WISE at 12 jum (blue 
source inside the 1 arcmin radius white circle). 

The emission in Herschel bands corresponding to 70 (blue), 
250 (green) and 500 jum (red) is shown as a color-composite im¬ 
age in Fig. [5] where several regions of interest can be seen. As 
can be noticed, the emission observed in the three different bands 
is similar, in the sense that it forms a shell-like structure, instead 
of having, as observed in several H ii regions, the emission at 70 
jum brighter in the interior of the bubble, where the hotter mate¬ 
rial is expected to be present, and the SPIRE emission, mapping 
the cooler dust, outside it (e .g. iRodon et alJl20ld : IPeharveng et alJ 
I 2 OIOI : [Anderson et al.ll2012h . In this context, the observed large- 
scale distribution of the infrared emission at the periphery of the 
ionized gas represents matter accumulated during the expansion of 
the ionized gas. Additional evidence supporting this conclusion can 
be seen by comparing the emission distribution among the differ¬ 
ent infrared bands, as shown in Fig [6| In this figure the different 
lines represent the mean value of the corresponding emission along 
concentric rings spaced by 1 arcmin and centred at (l,b) = (68H7, 
1?035). We divided the structure in two halves, the West and East 
sides, shown in the top and bottom panels of Fig. [b] respectively. 
The normalized emission was computed by dividing every point by 
the corresponding maximum value, 174.5,215.2 and 39.4 MJy/sr 
at 70, 250, and 500 jum respectively, for the East side; and 114.9, 
223.1, and 40.2 MJy/sr for the West half. From this comparison it 
can be seen that the emissions in the three bands present a simi¬ 
lar behavior and peak together at the location of the shell structure. 
It can also be noticed that in the West side the emission at 70 jum 
peaks twice. In this case the maximum located closer to the center 
of the structure is caused by the strong emission observed bordering 
the molecular cloud labeled as Cloud 1 in Fig.[T] 

The fact that the emission distributions show a similar struc¬ 
ture at 70 jum and at longer wavelengths suc h as 500 jum is consis- 
tent with the evolutionary path suggested bv ISmith & Contil (l2008h 
for a WNH star, where, as mentioned in the Introduction of this 
paper, the energy injected by the WR star impinges directly onto 
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the bubble created during the 0-type phase, making it possible that 
all the dust had been sweep-up into a shell structure. On the other 
hand, the spatial correlation between the infrared and CO emis¬ 
sion is noticeable, suggesting that the infrared emission originated 
in dense gas, where several CO, ^^CO, CH-i-, and Nil transition 
lines can take place and be responsible for most of the observed 
emission, as was d etected in the area of the Hii region Sh2-104 
(iRodon et al.ll20ld) . As can be seen in Fig. [3] the dense structure 
is a fragmented shell of collected material (e.g. cloud 3) with few 
probably pre-existing clumps (e.g. clouds 1 and 4) that clearly dis¬ 
torted the ionization front. 

Given that the physical conditions of the structure can be de¬ 
termined from the dust thermal emission, we have selected several 
regions (see Fig. [5]) for which we obtained the cold dust tempera¬ 
ture using the estimated Herschel fluxes, given in Table [2| To get 
the flux density in each band, the background, chosen far from the 
structure, was first subtracted. The regions were selected in order 
to sample dust related to different parts of the structure, i.e., the in¬ 
terior, the shell, the photo-dissociation region, the CHII region can¬ 
didate and cYSOs (see Section[5]), and the Hii region G68.14-^0.92 
located at a much larger distance. 

The spectral energy distributions (SEDs) obtained for the 12 
regions, normalized to the 160 jjm flux value, are shown in Fig.|7] 
An inspection of the distributions shows that all regions but three 
have a similar SED and are therefore expected to contain dust with 
similar properties. The regions that present a different SED are indi¬ 
cated with different colors and correspond to region 4 (green line), 
region 5 (red line) and region 7 (orange line). Erom Eigure[5]it can 
be seen that region 4 is associated with G68.14-^0.92 (the farthest 
Hii region, which is clearly visible in the infrared), region 5 is lo¬ 
cated in the interior of the shell and region 7 corresponds to the 
CHII region candidate (see Section[5]). 

To characterize the dust located in each region, the spectral 
energy distributions (SEDs) were fitted using the models created 
by DUSTEM, and a full description of the assu mptions and tests 
perfor med can be found in the original paper by ICompiegne et alJ 
(1201 ih . In a nutshell, DUSTEM creates spectral energy distribu¬ 
tions of the emission from interstellar dust covering a wavelength 
range from the mid to far-IR (approximately 3 to 1000 jim). The 
emitting dust is the result of three main components: polycyclic 
aromatic hydrocarbons (PAHs), amorphous carbon and amorphous 
silicates. The abundances and size distribution of each of these 
components and their interaction with the interstellar radiation field 
(ISRE), determine essentially the shape of the SED. Roughly, PAHs 
are expected to dominate the spectra between 3-15 jim, small 
amorphous carbon grains the 20-40 fim range and larger grains 
(amorphous carbon and silicates) the longer wavelengths (see e.g. 
ICompiegne et alJEoilL Fig.2). The present version of the code, al¬ 
though it is able to deal with ISO and Spitzer mid-IR data, needs 
to be modified to include WISE data. Therefore, in this study we 
are concentrating on the longer wavelength range covered by the 
Herschel Space mission. In terms of determining a dust tempera¬ 
ture representative of the medium, the large grains are the ones that 
can reach a thermal equilibrium with the ISRE, while smaller dust 
particles are stochastically heated, and therefore, a dust equilibrium 
temperature is meaningless. 

We thus estimated the dust properties of four different regions, 
the ones mentioned above which present a different SED (see Eig. 
[7]) and region 11 as a reference for the rest of the regions. Eor each 
region t he models provide the ISRE need ed to reach the thermal 
balance dMathis. Mezger & Panagial [l983h . and the temperatures 
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Figure 6. Averaged 70 //m (blue), 250 //m (green), and 500 jum (red) in¬ 
frared emission distributions as a function of angular distance from the cen¬ 
ter of the structure, (/, b) = (68? 17, 1?035). Top panel: East side of the struc¬ 
ture. Bottom panel: West side. 


of the three dust components (small carbon, large carbon and sili¬ 
cates). All the estimates are shown in Table [3] 

As can be seen in Table [3] the estimated temperatures for the 
Hii region G68.14-^0.92 and the candidate CHII (see Section|5]) are 
typical for H ii regions and higher than the ones found for the shell 
structure and its interior. On the other hand, the interior region is 
the coldest, with a temperature around 2 0 K, which is chyacter istic 
of the more diffuse interstellar medium dBoulanger et al.lll996h and 
consistent with the lack of bubble material in this region. As for 
the ISRE factor, it can be noted that it is also larger in regions 4 
and 7, as expected for regions near a strong radiation field, as is 
the case for these regions, since they host new stars (see Section|5])- 
On the contrary, for the interior of the shell (region 5) and the shell 
itself (region 11) the model yields a very low ISRE, suggesting that 
strong radiations fields are not involved in these regions. 


5 YOUNG STELLAR POPULATION 

As shown in the previous sections, it is clear that the interstel¬ 
lar medium surrounding WR 130 is strongly affected by the high 
UV radiation and mass loss rate of the star. As proposed by 
lElmegreenI dl998h . in this scenario it is highly possible that in¬ 
duced star formation is taking place. With this in mind, we look 
for primary tracers of stellar formation acti vity in the area makin g 
used of the IRAS Point Source Catalogue dBeichman et al. 198^ . 
the MSX Infrared Point Source Catalogue dEgan et~ l2003h .~the 
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Figure 5. Composite image of the region around WR 130. The image shows SPIRE 500 //m emission in red, SPIRE 250 jum in green and PACS 70 //m in blue. 
Outlined are the regions where flux densities were estimated. 


Table 2. Infrared flux densities estimated for regions 1 to 12. 


Region 

/(deg.) 

^(deg.) 


Flux Densities (Jy) 

22{]im) 70(jum) 160(jum) 

250(//m) 

35O0um) 

500(jum) 

1 

67.99 

1.03 

0.4 

0.7 

6.4 

32.5 

26.3 

13.4 

5.7 

2 

68.03 

1.05 

0.6 

0.9 

13.8 

20.6 

11.5 

4.4 

2.2 

3 

68.07 

1.07 

0.9 

2.2 

35.7 

55.3 

32.1 

16.0 

6.3 

4 

68.15 

0.92 

4.3 

16.1 

191.3 

212.1 

95.6 

36.2 

15.1 

5 

68.17 

1.04 

0.2 

0.3 

0.6 

5.8 

5.7 

2.9 

1.4 

6 

68.24 

1.12 

0.1 

0.2 

2.4 

6.6 

4.5 

2.2 

1.2 

7 

68.24 

0.96 

0.7 

5.5 

40.9 

23.5 

8.8 

3.3 

1.2 

8 

68.25 

1.18 

1.3 

1.8 

19.3 

86.0 

66.0 

30.2 

13.8 

9 

68.27 

1.12 

0.1 

0.3 

3.2 

8.2 

6.2 

3.1 

1.3 

10 

68.29 

0.85 

0.4 

0.8 

6.6 

32.7 

28.3 

13.8 

6.2 

11 

68.29 

0.98 

0.6 

1.1 

14.2 

38.2 

24.0 

10.5 

4.7 

12 

68.30 

1.07 

0.5 

1.2 

12.5 

25.5 

14.5 

6.4 

2.5 


WISE All-Sky Source Catalogue (IWright et aI.ll201Clh . and a pre¬ 
liminary source extraction from Herschel HiGAL data (Molinari et 
al. 2015, in preparation), using t he CUTEX source extraction and 
photometry dMolinari et al.ll201lh . Within a circular area of 20' ra¬ 
dius, centered at (/,/?) = (68°2, 1°0), a total of 29 IRAS, 10 MSX 
(having flux quality ^ > 1 in the four bands), 686 WISE (with 


signal-to-noise ratios greater than 7 and photometry flux uncertain¬ 
ties lower than 0.2 magnitudes, in the Wl, W2 and W3 bands), 
sources were found. 

To identify the young stellar object candidates (cYSOs) 
among the i nfrared sources, we adopted the classification scheni e 
described in ijunkes. Euerst & Reichl([l9^JLumsden et ahl (l2002l) . 

















8 S. Cichowolski et al. 


Table 3. Dust properties of the four selected regions. 


Region 

Ldust-small Carbons 

Tdust- large Carbons 

Ldust-silicates (K) 

ISRF 

4(068.14+0.92) 

34.4 

32.5 

22.3 

8.95 

5 (interior) 

20.1 

19.0 

13.7 

0.46 

7 (CHII) 

41.8 

39.4 

26.6 

25.5 

11 (shell) 

25.5 

24.1 

17.0 

1.72 



Figure 7. Infrared flux densities measured for the 12 regions shown in Fig. 
|5]at different wavelengths. The green, red, and orange lines correspond to 
regions 4, 5 and 7, respectively (see Fig.|5)- The black line was taken as the 
representative infrared energy distribution of the other 9 regions (shown in 
gray). 


and iKoenig et al.l (l2012h . for the IRAS, MSX and WISE data, re¬ 
spectively. For the Herschel data we have use d the prescription de¬ 
veloped by the HiGAL team (lEIia et alJl2013h , where essentially a 
source is considered a candidate YSO if it is detected at 70 fim and 
appears in at least 3 contiguous bands. 

In this way, we found that 2 IRAS sources, IRAS 19563-f3114 
and IRAS 19571-f3113, are cYSOs. As for th e MSX sources, ac¬ 
cording to the criteria of iLumsden et al.l (l2002h . we have identified 
massive young stellar object (MYSO) candidates and compact Hii 
region (CHii) candidates. In this region, among the 10 MSX listed 
sources we found 2 MYSO and 3 CHii candidates. We then ana- 
lyzed the listed sou rces from the WISE catalogue and, following 
iKoenig et al.l (l2012h . before attempting to identify the cYSOs, we 
first selected the non-stellar sources, such as PAH-emitting galax¬ 
ies, broad-line active galactic nuclei (AGNs), unresolved knots of 
shock emission, and PAH-emission features. A total of 119 sources 
were dropped from the list, and among the remaining 567 only 2 
were identified as Class I sources (e.g. sources where the IR emis¬ 
sion arises mainly from a dense infalling envelope) and 9 as Class 
II sources (e.g. sources where the IR emission is dominated by the 
presence of an optically thick circumstellar disk). In Fig.[8]we show 
the (W2 - W3) vs (W1 - W2) color-color diagram for all the un¬ 
contaminated sources, where the Class I and Class II sources are 
shown in red and green, respectively. Given that protostellar ob¬ 
jects with intermediate/high masses can be identified among the 
Clas s I sources by addit ionally requiring that the magnitude W3 
< 5 (iHiguchi et al .1120131) . and that neither of the detected 2 Class 



Figure 8. WISE color-color diagram showing the areas wh ere the Class 
I and Class II YSO candidates are located according to the IKoenig et al.l 
<20lA criteria. 


I sources satisfies this criterion, it seems that both are low mass 
protostars. 

In the far-infrared, using the Herschel data, one is systemat¬ 
ically detecting more deeply embedded objects; although is clear 
that those less embedded are likely to be detected at shorter wave¬ 
lengths as well. Unlike low mass protostars where there is a general 
agreement on their classification as function of their SED and evo¬ 
lutionary stage, from Class 0 to Class III, indicating t heir age and 
the c learance of their proto-stellar envelope (see e.g. lEvans et aP 
l2009h . for massive YSOs is not the case. It is possible, however, 
und er the ’mass accr etion paradigm’ of massive star formation (see 
e.g. iTan et al.ll20l4 and references therein), to develop a similar 
scheme where the youngest mor e embedded massive YSOs are 
identified as Class 0 and so on (iMolinari et al.l |20()3) . Observa- 
tionally, the HiGAL team has developed a method to select YSO 
’candidates’ that allows such classification and places the sources 
in a bolometric luminosity versus mass envelope diagram, provid¬ 
ing an idea on their evolutionary stage. The scheme, as mentioned 
above, requires a detection at 70 fim and three contiguou s bands 
(to create a SED); for further details see (lElia et al.ll2013h . In this 
study we have used such scheme, but we have gone a step fur¬ 
ther by trying to fit a SED model to data when possible, and us- 
ing the Robitaille-Whitney H yperion grid of models for such task 
jRobitaille & Whitnevll2014h . 

The location of all the cYSOs found are indicated in Fig. [9] 
and their infrared properties listed in Tabled From Fig. [9] it is ev¬ 
ident that star formation is taking place in the local environment 
of WR130, and that the cYSOs are found preferentially onto the 
molecular clouds labeled as Cloud 1, 2, 3, and 4, in Fig. [T] It is 
clear that the spatial distribution of the cYSOs follows the distribu¬ 
tion of the high column density matter. 
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Table 4. IRAS, MSX, WISE and Herschel sources found to be YSO candidates. 


IRAS sources 


# 

Designation 

(l,b) 

Fi 2 [Jy] (Qi 2 ) 

F25[Jy] (Q 25 ) 

F6o[Jy] (Q 60 ) 

Fioo[Jy] (Qloo) 

Notes 

1 

19563+3114 

68.07, 1.07 

1.690 (3) 

4.411 (3) 

56.06 (3) 

105.2 (2) 

Cloud 1 

2 

19571+3113 

68.146, 0.92 

3.682 (1: upper limit) 

15.050 (2) 

173.80 (3) 

399.9 (2) 

G68.14+0.92 

MSX sources 

# 

Designation 

{l,h) 

FgiJy] (Qg) 

Fi 2 [Jy] (Q 12 ) 

Fi4[Jy] (Qi4) 

F 2 i[Jy] (Q 21 ) 

Notes 

3 

G068.0331+01.0543 

68.03, 1.05 

0.74298 (4) 

1.3394(1) 

0.85034 (1) 

2.7624 (1) 

CHII, Cloud 2 

4 

G068.1492+00.9247 

68.149, 0.92 

1.08850 (4) 

1.9132 (3) 

1.42180 (3) 

3.0316(3) 

CHII, G68.14+0.92 

5 

G068.2383+00.9596 

68.238, 0.9596 

0.96607(4) 

1.4840 (3) 

1.47640 (3) 

8.0493 (4) 

CHII, CO also at -62 km s“^ 

6 

G068.3277+01.0201 

68.32+01.02 

0.60381 (4) 

2.30430 (4) 

2.8046 (4) 

4.9903 (4) 

MYSO, Cloud 3 

7 

G068.1424+00.9197 

68.14, 0.919 

0.37255 (4) 

0.69779 (1) 

1.2673 (3) 

3.4386 (3) 

MYSO, G68.14+0.92 





WISE sources 




# 

Designation 

(l,b) 

W1 [mag] 

W2 [mag] 

W3 [mag] 

W4 [mag] 

Notes 

8 

J195819.33+313643.5 

68.266, 1.20 

14.366 

12.894 

8.784 

7.292 

Class I, Cloud 4, CO also at -60 kms“^ 

9 

J195839.60+313404.4 

68.266, 1.11 

13.436 

12.118 

7.874 

5.252 

Class I, Cloud 4 

10 

J195824.18+313643.7 

68.27, 1.185 

9.115 

8.575 

6.160 

3.996 

Class II, Cloud 4 

11 

J195833.72+313453.6 

68.267, 1.140 

11.840 

11.493 

7.501 

5.581 

Class II, Cloud 4 

12 

J195820.08+312011.3 

68.033, 1.053 

9.388 

8.388 

4.298 

2.370 

Class II, Cloud 2 

13 

J195925.21+312333.1 

68.20, 0.88 

14.019 

13.607 

8.751 

7.244 

Class II, no CO 

14 

J195832.30+313304.2 

68.238, 1.129 

13.383 

12.727 

8.000 

5.588 

Class II, Cloud 4 

15 

J195854.39+311831.1 

68.073, 0.935 

10.025 

9.456 

6.587 

4.210 

Class II, no CO 





Herschel sources 



# 

Designation 

(l,b) 

Fvo [Jy] 

F 160 [Jy] 

F 250 [Jy] 

F 350 [Jy] 

F 500 [Jy] 

16 

238687 

67.973, 1.022 

0.329 

0.761 

1.611 

- 

- 

17 

238749 

68.008, 1.049 

0.785 

1.689 

5.746 

4.558 

2.004 

18 

238786 

68.026, 1.062 

1.3 

2.864 

6.511 

12.94 

2.957 

19 

238789 

68.028, 0.983 

1.276 

5.812 

5.762 

3.265 

1.663 

20 

238797 

68.033, 1.054 

11.424 

16.491 

12.655 

2.347 

2.293 

21 

238835 

68.057, 1.076 

0.455 

3.849 

6.204 

- 

- 

22 

238863 

68.069, 1.062 

15.427 

20.963 

13.061 

10.912 

6.197 

23 

238880 

68.078, 1.074 

74.122 

49.372 

36.245 

18.919 

9.604 

24 

238886 

68.08, 0.959 

0.804 

3.267 

2.641 

1.705 

2.219 

25 

238887 

68.08, 0.966 

0.891 

2.022 

4.1 

- 

- 

26 

238888 

68.080, 1.06 

1.053 

3.049 

1.564 

- 

- 

27 

238953 

68.124, 0.939 

1.478 

4.173 

3.185 

5.83 

2.375 

28 

238963 

68.130, 0.921 

0.551 

1.354 

2.001 

- 

- 

29 

238979 

68.138, 0.933 

2.115 

3.978 

1.344 

- 

- 

30 

238993 

68.143, 0.897 

8.92 

14.102 

13.491 

0.501 

- 

31 

239001 

68.145, 0.905 

1.717 

10.201 

6.913 

12.91 

6.887 

32 

239003 

68.146, 0.922 

5.407 

27.825 

- 

- 

- 

33 

239007 

68.148, 0.897 

1.284 

3.301 

- 

- 

- 

34 

239008 

68.148, 0.915 

3.876 

3.983 

16.76 

- 

- 

35 

239012 

68.150, 0.928 

16.703 

70.192 

35.56 

- 

- 

36 

239013 

68.150, 0.922 

8.439 

46.275 

35.76 

21.896 

11.192 

37 

239094 

68.195, 1.125 

2.781 

13.222 

9.275 

4.634 

2.413 

38 

239192 

68.24, 0.959 

37.411 

18.088 

15.479 

5.295 

2.463 

39 

239252 

68.266, 1.116 

1.249 

4.231 

2.361 

- 

- 

40 

239271 

68.276, 1.185 

0.955 

1.444 

1.758 

0.639 

- 

41 

239279 

68.279, 0.948 

0.739 

6.536 

12.99 

8.077 

4.812 

42 

239302 

68.294, 1.065 

0.73 

1.843 

6.333 

1.531 

1.218 

43 

239366 

68.323, 1.076 

1.567 

13.42 

7.726 

11.75 

4.159 

44 

239370 

68.327, 1.075 

5.431 

11.735 

16.903 

- 

- 

45 

239371 

68.328, 1.02 

4.212 

0.872 

0.808 

- 

- 
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Figure 9. Herschel image at 70 jim with the location of the cYSOs superim¬ 
posed. Green, blue, and red asterisks correspond to IRAS YSO, MYSO, and 
CHii candidates, respectively, while black and red circles indicate the loca¬ 
tion of WISE Class I and Class II candidate sources, respectively. Herschel 
sources are indicated by red crosses. 


As can be noticed, the IRAS source IRAS 19571+3113 (# 
2 in Table |4]), two MSX sources (#4 and #7), and 9 Herschel 
sources (#28 to #36) are l ocated onto the H ii region G68.14+0.92, 
which, as pointed out bv Cichowolski et al. ( 200ll), is no t related 
to WR130. Moreover, iBronfman. Nvman & Mavl (Il996h in their 
survey of the CS(2-1) emission toward IRAS point sources in the 
galactic plane, observed, in direction to IRAS 19571+3113, CS(2- 
1) emission at -62.6 kms“^ suggesting that a high density molec¬ 
ular clump is associated with G68.14+0.92 and may harbor new 
stars. Given that this region is far away from the area under study, 
these 12 infrared sources will not be considered further. 

As regards to the MSX source candidate to be a CHII 
G068.2383+00.9596 (# 5 in Table lU, it coincides with the Her¬ 
schel source #38 and, from Fig.|3 it can be seen that, in projection, 
their location is very close to WR 130. This source is observed 
both at radio and infrared wavele ngths and it is includ ed in the 
Red MSX Source (RMS) survey (lUrauhart et al.ll2008h . Accord¬ 
ing to these authors, G068.2383+00.9596 is associated with the 
infrared source IRAS 19572+3119, several 2MASS point sources, 
the WISE source 5087 (with estimated flux values only at 12 and 
22 jum ), and the radio source G068.2389+00.9592 (lUrauhart et al.l 
l2009h . for whic h a flux intensity of 6 .1 mJy is given at 5 GHz. On 
the other hand, lUrquhart et al.l J201lh detect neither ammonia nor 
water maser emission in the direction to this infrared source, and 
cannot thus suggest any kinematical distance to it. An inspection of 
the CO data cube shows, however, that CO emission in this direc¬ 
tion is only detected at about -62 km s“\ as shown in Fig.[T0l sug¬ 
gesting a further distance for the CHII candidate, of about 12 kpc, 
similar to the distance of the Hii region G68.14+0.92, indicating 
that we have located an active star formation region that is located 
far beyond the solar circle. In this case, the infrared sources #5 and 
#38 can not be associated with WR 130. On the other hand, an in¬ 
spection of the entire CO data cube reveals that sources #13 and 
#15 (see Table lU are not immersed in any molecular cloud, indi¬ 


Glon: 68.2375° 



Figure 10. CO emission in direction to the MSX source 
G068.2383+00.9596. 


eating that, if they are YSOs, either they have already moved away 
from the place where they were formed or that they dissipated it. 

From the remaining cYSOS appearing in Table ID it is impor¬ 
tant to mention that CO emission in their direction is only observed 
at the radial velocity interval between -9 and -15 kms"^ coinci¬ 
dent with the radial velocity of the molecular structure, suggesting 
an association between them and the molecular clouds. The excep¬ 
tions are Source #8, for which CO emission is also detected at -60 
km s~^ and Sources #22 and #26, in which direction there is CO 
emission also at-40 km s“^ 

Among the cYSOs, two groups of sources are of particular 
interest. As can be seen from Fig. |9l four sources, #1-23 (since 
sources #1 and #23 spatially coincide), #21, #22, and #26 are seen 
projected on the periphery of the CO cloud named as Cloud 1 in 
Fig. [T] while three sources, #3-12-20, #17 and #18, are locatated 
onto Cloud 2. In both cases the cYSOs are located on the side nearly 
facing the star WR 1 30. This suggests that they could be bright- 
rimmed clouds (e.g., |Morgan. Urquhart & Thornpson|[2009h being 
excited by the ionizing radiation from WR 130, in whic h case one 
could expect to And a radio continuum counterpart fe.g. JSchwartzl 
Il985h . An inspection of Fig. HU shows indeed the presence of a 
weak radio continuum source possibly barely resolved (it is hard to 
clearly separate it from its surrounding emission) with the CGPS 
beam, thus suggesting a size of approximately I'or less. A peak 
1420 MHz flux density of ~ 10+1 mJy is obtained with an es¬ 
timated total flux density of 18 + 3 mJy in the border of Cloud 1 
and a peak flux value of about 2 mJy is estimated for the emission 
bordering Cloud 2. 

Assuming the radio cont inuum emission to be free-fr ee emis¬ 
sion, we can use the models of iMezger & HendersonI (Il967h to infer 
th e properti e s of th e ionized gas, following the procedure outlined 
byi ineaultl (Il998h . Since we do not know a priori the actual dis¬ 
tribution of the gas, we take the simplest possible model for the 
source, namely that of a sphere with constant density. 

Then the electron density He, mass of ionized gas M and tota l 
number Nu of UV ionizing photons per second fe.g. lChaissorJl976l) 
are given by: 
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Table 5. Estimated parameters for the boundary ionized layers related to clouds 1 and 2 . 


Cloud 

rig (cm ^) 

Mi (Mo) 

Nu (10‘*o s-‘) 

<R (10'*'^ s-1) 

Pilk (10'’ cm-3 K) 

PmlktXd' cm-^K) 

1 

29 

2.2 

2.0 

3.0 

1.9 

2.0 

2 

12 

1 

0.4 

3.5 

0.6 

6 



Galactic Longitude 


Figure 11. 1420 MHz emission distribution in the area of Sources #1-23 and 
#3-12-20 (both indicated by black asterisks). Red contours are at 9.2, 9.6 
and 10 K. Black contours correspond to the averaged CO emission between 
-2.35 and 4.25 kms“^ from 1 to 2 K with a step of 0.2 K. The location of 
WR 130 is indicated by the blue cross. Herschel cYSOs are indicating by 
their source numbers (see Table |4). 

M/Me = 0.386 «2 ^ 

Nu = 0.76 X 10'*" 

where Sy is the measured flux density in Jy, Oca the measured 
half-power width in minutes of arc, T 4 the electron temperature 
in units of 10"^ K, a{y, Tg) the Gaunt factor (of order unity), Dkpc 
the distance in kpc and vghz the frequency in GHz. The parame¬ 
ters ui and U 2 are of order unity for a spherical source model. As 
the observing frequency appears with a very low exponent, all fre- 
quen cy factors are also of or der unity. Assuming a distance of 5 
kpc jCichowolski et al]l200lh . a temperature of 10"^ K and using 
Ogsl ~ 1 (valid for a source unresolved with a beam of the order of 
1 '), the above three equations thus simplify to 

n, -2905®'^ cm-^ 

M/Mq ~ 22Sf, 

Nu ~ 1.9x 10 ‘**Sv s“‘. 

It should be emphasized at this point that the above esti¬ 
mate for He is likely an underestimate since the source geome¬ 
try is that of a thin layer of ionized gas on the star-facing side 
of the cloud. Since the electron number density is of the form 
n, oc e-y oc (volume) we expect rig to scale roughly as 77 
if we assume the thickness of the layer to be given by rj times the 
cloud radius. 


Using the measured peak flux density of 10 mJy and 2 mJy for 
the border of cloud 1 and 2 respectively, the above numbers reduce 
to the values given in Table|5] As both sources are essentially unre¬ 
solved at V the estimated densities and ionized masses are a lower 
and upper limit, respectively. Their corresponding maximum size 
is on the order of 1.5 pc for a distance of 5 kpc. 

Assumi ng that WR 130 is the star providing the ionizing radi¬ 
ation, which der Huchl (l200lh gives as Awr = 5.0 x 10"^^ s“^ 
we can calculate the number of ionizing photons intercepted 
by the radio continuum features bordering cloud 1 and 2. For this 
estimate we considered that the border of cloud 1 is located some 
12' or 17.4 pc from WR 130 (we do not make any correction to 
the distance or cloud size due to projection effects as these cancel 
out), and the border of cloud 2 is a little bit further from WR 130, 
at about 12'5 or 18 pc. We And that a fraction fl ~ 0.006 and 
/2 ~ 0.007 of the photons from WR 130 would be intercepted for 
each cloud, so that ~ 3.0 x 10"^^ s“^ and ~ 3.5 x 10"^^ s“^ . 
It would seem, at least at first sight, that WR 130 alone is sufficient 
to power the observed radio continuum emission. 

We thus have two molecular clouds, labeled Cloud 1 and 2, 
presenting a bright rim and each harboring cYSOs. This together 
with the fact that the observed cometary shapes of the bright rims 
are quite apparent with their tips pointed towards the location of 
WR130 (see Fig. [9]) suggest the possibility that the cYSOs were 
formed under the radiative driven implosion (RDI) mechanism 
(lBertoldilll989h . This mechanism takes place when a photoioniza¬ 
tion shock is driven into the molecular gas. To ascertain whether 
this is the case, it is necessary to compare the pressure of the ionized 
layer bordering the molecular cloud with the molecular pressure. If 
the cloud is under-pressured with respect to the ionized layer or the 
pressures are similar, the ionization front will be able to propagate 
into t he cloud and then modify its evolution (iLefloch & Lazarem 
11994 . 

The pressure in the ionized boundary layer can be evaluated 
from the electron density as P/ = 2 p/Cp where pt is the density 
in the boundary layer, pt - UgniH and C/ is the sound speed in 
the ionized gas (ci ~ 11.4 kms“^). As mentioned above, this pres¬ 
sure should be compared with the pressure of the molecular cloud, 
which can be estimated, assuming that the thermal component can 
be neglected, from the turbulent velocity dispersion, cr^, and the 
molecular density, p^, as Pm - cr^ Pm, where cr^ may be written 
as cr^ = Al >^/(8 In 2 ), Av being the observed velocity line width 
of the molecular cloud gas. The results are given in Table|5] Before 
comparing the obtained values it is important to mention that the er¬ 
rors involved in these estimates are large. As mentioned before, the 
electron densities obtained for the ionized gas are lower limits, and 
thus the estimates of P/ are also underestimates. On the other hand, 
since the molecular pressure depends on the geometry adopted to 
estimate the molecular density, the CO-to-H 2 conversion factor X, 
and the distance adopted for the structure, its value has at least a 
60% error. Bearing this in mind, from Table |5] we can see that the 
pressures P/ and Pm obtained for cloud 1 are similar while in the 
case of cloud 2 the molecular pressure is higher than the pressure 
of the ionized gas, suggesting that the shock could not yet affect 
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the molecular gas. However, given the large errors involved in the 
estimates, with the data at hand we cannot completely rule out the 
possibility that the shock front has already penetrated into cloud 2. 

The second group of sources of interest in Fig. |9] is located 
onto Cloud 4 (see Fig.[T]). This group consists of five WISE sources 
(two of them coincide with Herschel sources), two classified as 
Class I YSOs and three as more evolved. Class II, objects (see 
Table |4]). An inspection of Fig. |4] clearly shows, more strikingly 
at 250 yum, two jet-like features spatially coincident with the bor¬ 
der of molecular Cloud 4. Moreover, two of the cYSO, identified 
as sources #9, #39 and #14 in Table IH are located onto the tip 
of these structures. A closeup of this region is shown in Fig. [12] 
The cometary shape of the infrared structures suggests that the RDI 
mechanism could be at work in this molecular cloud, but the low 
angular resolution of both the molecular and radio continuum data 
prevents carrying out the analysis in this region. 

Among the rest of the cYSOs presented in Table IH it can be 
seen that four sources (#6-45, #42, #43, and #44) are located onto 
Cloud 3 (see Fig.[T]), while the rest are distributed onto the general 
CO shell structure. 

5.1 Physical Properties of cYSOs 

To characterize the cYSOs identified located in projection onto the 
molecular clouds, we constructed their SEP using th e grid of mod¬ 
els and fitting tools of iRobitaille et Hl bOOfiLlaM . The SED fit- 
ting tool fits the data allowing the distance and external foreground 
extinction as free parameters. We gave a distance range of 4.5 - 5.5 
kpc. The visual absorption range for each source was set making 
use of the Hi and CO data cubes to estimate their column den¬ 
sities at the distance of th e structure, using the relation N (HI) -h 
2 N(H 2 ) ~ 1.9 X 10^^Ay teohlin. Savage & Drake! 1 197^ . where 
N(HI) is the Hi column density and is given by N(HI)(cm"^) = 
1.824 X 10^^ r(K) Av (km s"^), and N(H 2 ) is the molecular column 
density, given by N(H 2 )(cm"^) = 2.3 x 10^*^ T (K) Av (kms"^), 
where, in both cases, T is the brightness temperature and Av is 
the velocity interval between the structure and the observer. The 
adopted range for each source is given in the last column of Table 

m 

Eor each cYSO besides using the photometry given in Table 
|4]we search their counterpart in other infrared ranges, to use in the 
fitting as many points as possible to better constrain the estimated 
parameters. Thus, for sources #3-12-20, #11, #14 and #15, JHK 
photometry was obtained from the 2MASS catalogue. Eor source 
#1-23 the flux given in the MSX catalog at Band A was also con¬ 
sidered. Eor source #3-12-20 out of the four MSX bands observed 
only Band A was considered due to the bad quality of the other 
three, and the IRAS fluxes obtained from the IRAS catalogue were 
used as upper limits. 

Eor each source, based on all the fitted models satisfying 
- xlcst < 2 A, where xl^st goodness-of-fit parameter for 

the best-fit model and N is the number of input observational data 
points, we estimated some parameters such as the stellar mass, stel¬ 
lar temperature, stellar age, envelope accretion rate, envelope mass 
and total source luminosity. To estimate the parameters, we care¬ 
fully inspect every distribution considering all the models that sat¬ 
isfy the criteria mentioned above. The results are listed in Table [6] 
We note that, with the data at hand, even when the data points can 
be fitted with great goodness, just a few parameters are well con¬ 
strained and in some cases only lower and/or upper limit values can 
be given. 

As can be seen in Table [6] all the sources but five are in a very 



^ ^ ^^ ^^^^^^^^^^^ 

68.35° 68.30° 68.25° 68.20° 68.15° 

Galactic Longitude 


Figure 12. 250 jum emission distribution in the area of the molecular cloud 
4. The five WISE sources located onto this cloud are indicated with aster¬ 
isks and with their corresponding numbers (see Table|4). Gray contours are 
from 0.4 to 1.6 K with a step of 0.4 K and correspond to the averaged CO 
emission between -7.18 and -0.7 kms“k The location of WR130 is indi¬ 
cated by the blue asterisk. 


early stage of their formation, with ages lower than 10^ yr, still hav¬ 
ing high envelope masses, of the order of hundreds, and envelope 
accretion rates in the order of 10“^ - 10"^ M© yr“^ Similarly, most 
of the cYSOs have low/intermediate masses, lower than 7 M©, al¬ 
though these values will probably increase considerably since they 
are still accreting lot of mass from their envelopes. 

As mentioned before, Eig.j^ shows that sources #1-23, #21, 
#22, and #26 are located within the molecular cloud 1, while 
sources #3-12-20, #17 and #18 are within cloud 2. An inspection 
of the ages estimated for these cYSOs shows that in cloud 1 the 
age gradient is as expected for sources #1-23, #22 and #26, if the 
RDI mechanism took place in this cloud, in the sense that the oldest 
source (#26, among these three) is located closest (in projection) to 
the WR star, then is source #22 and then #1-23. The problem is with 
source #21, which, since it is located further from WR 130 should 
be younger. A similar situation occurs for the three sources located 
within cloud 2, where the ages of sources #17 and #18 are larger 
than the one of #3-12-20. Looking at the SEDs fitting we found 
that for these three sources particulary the fitting is not satisfactory. 
More data at other wavelengths would be useful to better estimate 
the age of these sources. 

Three of the sources that have higher ages are sources #8, #10- 
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Table 6. SED parameters obtained for the cYSOs. 


# 

Stellar Mass 

Mo 

Stellar Temperature 

K 

Stellar Age 
10^ yr 

Envelope mass 
Mo 

Envelope Accretion Rate 
Moyr-l 

Total Luminosity 
105 Lo 

Av 

mag 

1-23 

8 

4130 

1.3 

476 

5.6x10-* 

2.5 

4-10 

3-12-20 

4 

4143 

1.7 

20 

0.002 

0.3 

4-8 

6-45 

6.8 

4435 

21 

0.4 

3.7 xlO-5 

0.9 

2-5 

8 

3-6 

11000- 18700 

> 1200 

< 7 X 10"^ 

0 

0.05-0.97 

2-4 

9-39 

1.2 

3967 + 15 

1.1-1.3 

69-72 

(0.4±0.1)xl0-* 

0.04-0.05 

2-6 

10-40 

4.6 

11270 

836 

10.6 

1.5 xlO"'' 

0.4 

3-7 

11 

5.3+ 0.4 

4688 +118 

95-185 

8-13 

(0.3- 1.0) xlO-'* 

0.1-0.2 

2-6 

13 

4.2 ±0.3 

14855 +600 

2200-4100 

< 9 X 10-^ 

0 

0.26 +0.06 

2-6 

14 

1.9 

4320 

112 

12 

9 xl0“^ 

36 

2-7 

15 

3.5 

11820 

1954 

10"^ 

0 

0.2 

2-4 

16 

1.3+ 0.3 

4092 + 94 

7- 14 

97 + 17 

(50 ±7) X 10-5 

0.03 + 0.01 

5-8 

17 

3.4 

4390 

39 

251 

10-5 

0.09 

6-9 

18 

5.5 

4433 

28 

760 

2 xlO-5 

0.2 

7-10 

19 

5.5 

4433 

28 

760 

2 xlO-5 

0.2 

4-7 

21 

2.6 

4200 

7.6 

274 

10-5 

0.1 

6-9 

22 

2.3-2.5 

4215 + 22 

8-15 

223 - 274 

10-5 

0-07-0.10 

6-9 

24 

2.3-3.4 

4200 - 4393 

6.4-38.7 

223 - 345 

10-5 

71 - 112 

4-7 

25 

1.6-3.4 

4114-4393 

6-39 

165-251 

10-5 

0.05 - 0-09 

4-7 

26 

2.5 ± 0.8 

4254 ± 125 

23 ±9 

14 - 262 

(7±2)xl0-* 

0.08 ± 0.02 

4-8 

27 

5.5 

4433 

28 

760 

2 xlO-5 

0.2 

3-5 

37 

4.9 

4380 

6-20 

540 - 840 

2 xlO-5 

0.25 - 0.30 

4-6 

41 

4.6 

4208 

4.7 

1241 

2 xlO-5 

0.3 

5-8 

42 

2.3-3.4 

4200 - 4393 

7.6-38.8 

223 - 274 

10-5 

0.07-0.10 

2-5 

43 

4.6 

4208 

4.7 

1241 

2 xlO-5 

0.3 

5-8 

44 

4.6-8.1 

4423 + 189 

<30 

1142 + 228 

2 xlO-5 

0.5 ± 0.2 

5-8 


40 and #11, all lying in projection onto Cloud 4. In this case, if the 
RDI mechanism is responsible for the star formation in the cloud, 
given the location of these three cYSOs with respect to WR 130, 
they would have the lowest ages. Since this is not what we obtained, 
we suggest that either another mechanism is acting in cloud 4 or 
that both sources are not related to this cloud. As mentioned before, 
CO emission is also detected at -60 km s~^ in direction to source 
#8. Also sources #13 and #15 are among the oldest cYSOs, which 
is consistent with the fact that, as mentioned before, they are not 
inmersed in any molecular gas. 

The simultaneous presence of envelope-dominated YSOs 
(typical age around 10"^ years) and a WR star in the same region 
suggests that different episodes of star formation have occurred in 
the region. Since WR 130 is an evolved star, we can speculate that 
the formation of young sources in the region might have been trig¬ 
gered by the expanding bubble. 


6 SUMMARY 

The previous observations by ICichowolski et alJ (l200lh and those 
reported here allow us to describe a consistent picture of the 
recent evolution of WR 130 and of its effect on its immediate 
environ ment on large as well a s small scales. The intitial re¬ 
sults of ICichowolski et all (l200lh showed that the continuum ra¬ 
dio emission and H i kinematical analysis were in total agreement 
with the generally accepted picture of stellar wind e volution (e.g. 
ICastor. McCray & Weaveiil 19751 : 1 Weaver et al.lll977h . In this paper 
we have gone one step further and used recent high-quality IR and 
CO data to study for the first time the molecular and dust properties 
of the ISM local to this Wolf-Rayet star. 

The timescales involved are consistent with the wind of the O 


star phase having formed an initial cavity in the surrounding ISM. 
The wind of the star now in its WR phase seems to be interacting 
directly with the cavity wall as there is no sign of massive ejec¬ 
tion having taken place in the recent history of the star. If any such 
ejection took place, it is likely to have been a minor event. In all 
probability WR 130 is now in a pre-LBV phase. 

On large scales, an inspection of the CO data cube reveals 
the presence of several molecular clouds which have a good mor¬ 
phological agreement with the ionized gas. The CO structures are 
observed within the velocity range between -11 and 3 km s“^ com¬ 
patible w ith the velocity of the Hi structure and ionized gas as ob¬ 
tained by ICichowolski et ^ (l200lh . Among the observed molecu¬ 
lar clouds, four have been identified and characterized, over which 
several cYSOs are located. 

At infrared wavelengths, the large-scale distribution of the 
WISE emission at 12 /dm clearly shows the presence of a PDR bor¬ 
dering the HII region. In Herschel bands a shell like structure is no¬ 
ticeable, showing a good morphological coincidence with the CO 
structures, suggesting that the dust coexists with the molecular gas. 
This is consistent with dust temperatures of about 24 K obtained 
for the shell structure. 

On smaller scales, we see clear evidence of gas and dust being 
in the process of forming new stars, together with an accompanying 
erosion of the associated cloud, as a result of the combined action 
of the stellar UV radiation and wind of WR 130. This appears most 
clearly on the high-resolution WISE images which show morpho¬ 
logical features (e.g., oriented pillars or thin filaments) generally 
seen around young massive stars or stellar clusters and which are 
indicative of stellar formation. 

The presence of compact sources at the tip of some of these 
structures identifies these as probable YSOs. Two groups of such 
objects have been identified. In the first case, two of the compact 
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IR sources have an associated weak radio source located on the 
star-facing side of their associated CO cloud. Using the available 
information, we have shown that these two cases are consistent, 
one formally and the second one marginally, with a formation via 
the RDI mechanism. The other sources, which have a more elon¬ 
gated or pillar-like appearance, do not have associated radio con¬ 
tinuum counterparts and a similar analysis cannot be carried out. 
A different mechanism involving the curvature of a dense expand¬ 
ing shell leadin g to the formation of st a ble pill ar-like structures has 
been proposed dTremblin et al.ll2012al lR l201^ and could very well 
be a promising avenue to explain the observed features. However 
such an analysis is clearly out of the scope of this paper. 

Based on different colour criteria, an analysis of the YSO can¬ 
didates in the vicinity of WR130 has allowed us to estimate some 
of their physical parameters and revealed that most are in a very 
early phase of their formation, still accreting mass, and that they 
are mainly low mass stars. An estimate of their age shows that, 
whereas the age gradient is as expected in some cases, there are 
nevertheless cases where the situation is uncertain. It is not clear 
whether this effect is real or the result of a poor fitting of the SEDs. 

In summary, the morphological correlation between ionized 
and neutral gas as well as the ring-like dust emission confirm the 
high degree of interaction between the wind and UV photons from 
the WR star and the surrounding molecular clouds. In addition, the 
presence of numerous cYSOs within the molecular gas suggests 
that an active star formation process is taking place in the environs 
of this star. 
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